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Abstract 
(N 

>, 
v<0 . Standard Model singlet neutrinos propagating in extra dimensions induce small Dirac 

"^ I neutrino masses. While it seems rather unlikely that their Kaluza-Klein excitations 

f^^ ■ directly participate in the observed neutrino oscillations, their virtual exchange may 

O '' lead to detectable signatures in future neutrino experiments and in rare charged lepton 

processes. We show how these effects can be described by specific dimension-six effective 

operators and discuss their experimental signals. 

.^ 1 Introduction 

X. 

H , The hypothesis that Standard Model (SM) singlet fields propagate in extra dimensions leads to striking 
results. When applied to the graviton, it allows to lower the quantum gravity scale down to few TeV [1, 2], 
suggesting a new scenario for addressing the Higgs mass hierarchy problem. It is also natural to consider 
the case of "right-handed neutrinos" {i.e., fermions without SM gauge interactions) propagating in extra 
dimensions. The smallness of the neutrino masses, of the Dirac type, could in fact be a manifestation of 
this hypothesis [3, 4, 5, 6]. 

If the radius of the compactified dimensions is very large, R ^ eV~ , Kaluza-Klein (KK) modes of 
right-handed neutrinos would significantly participate in neutrino oscillations. However, KK interpreta- 
tions of the atmospheric and solar neutrino puzzles are disfavoured by the following arguments: 

• A KK tower of sterile neutrinos gives rise to active/sterile oscillations at a small Am^ ~ l/R? 
only in the case of a single large extra dimension. In the case of more large extra dimensions, the 
active/sterile mixing is not dominated by the lightest KK modes, and the mixing with the heaviest 
KK modes does not lead to oscillations. One could argue that right-handed neutrinos in a single 
very large extra dimension can be effectively obtained if one of the extra dimensions happens to 
be much larger than the other (s). However, the fact that infrared effects are dominant in this case 
destabilises the hierarchy [7]: the Newton or Coulomb potential grows with the size R of the largest 
dimension. 



'On leave from Dipartimento di Fisica deU'Universita di Pisa and INFN, Italy. 



• Even if one forgets about the hierarchy problem, there are severe bounds from supernovse obser- 
vations [8]. One needs to prevent resonant neutrino conversion in supernovae by choosing a smah 
radius R<,1/ MeV or by adding an ad-hoc 5-dimensional mass term m ^ MeV for the right-handed 
neutrino(s) responsible for "atmospheric" oscillations. The latter still allows to build models for 
the solar anomaly roughly compatible with supernovae bounds [9]. 

• Finally, the recent results from the SuperKamiokande [10, 11] and SNO [12] experiments provide 
strong indications for z^,- appearance in atmospheric oscillations and for u^^Vr appearance in solar 
oscillations. In both cases, a sterile interpretation is now strongly disfavoured by experiments. 

In this paper we explore the phenomenology of more promising models with 5 > 2 extra dimensions 
and radii which are not larger than what is required to reproduce the gauge/gravitation hierarchy. The 
neutrino puzzles are solved by "normal" (active) oscillations, but the presence of the heaviest KK neu- 
trinos can still lead to small but detectable effects in neutrino flavour transitions. After integrating out 
the heavy KK modes, we obtain an effective Lagrangian that contains massive Dirac neutrino states and 
a specific set of non-renormalizable operators. Since the dominant effects come from the heaviest KK 
states, the coefficients of these operators can only be estimated by introducing an arbitrary ultraviolet 
cut-off. 

At tree level one only obtains the dimension-six operators 

Aree = eij 2V2Gf{H^ U)iB{HLj), (1.1) 

where Li are the lepton left-handed doublets, H is the Higgs doublet, and eij = e*j are dimensionless 
couplings. The presence of these effective operators leads, for example, to potentially large flavour 
transitions P{ui — > Vj) ~ |e?| at 0{L^) {i.e., at very short baselines L <C Ey/ /^m?'), and CP-violating 
effects at 0{L}) (rather than at 0{L'^) and 0{L^) as in ordinary oscillations). Since this peculiar tree 
level operator only affects neutrinos, potentially detectable effects, especially at a neutrino factory, are 
not already excluded by bounds from rare charged leptons processes, like r — > fj,j, fi — > 67, fi — > eee, 
etc. We will show, however, that Eq. (1.1) cannot explain the LSND anomaly [13], due to the present 
constraints from searches for ^ — > 67. 

At one-loop level, other operators are generated, giving rise to rare muon and tau processes that 
violate lepton flavour. Again, the coefficients of these operators are cut-off dependent and can only be 
estimated. However, in minimal models, their flavour structure is directly related to the physical neutrino 
mass matrix, giving predictions for rare muon processes in terms of neutrino oscillation parameters. This 
is in contrast with other cases of physics beyond the SM. In supersymmetry, for instance, the rates for rare 
muon processes are perturbatively calculable, but their relations with neutrino oscillations parameters 
are strongly model-dependent and can vary by many orders of magnitude. 

Many of the effects studied here have already been considered in previous analyses [14, 15] as due to 
mixing between ordinary neutrinos and the whole tower of KK states. The equivalent language of effective 
operators we employ allows the study of different models {e.g. large and warped extra dimensions) in 
a unified framework, the discrimination of what is really computable from what can only be estimated, 
and a more transparent identification of all potentially interesting experimental signatures. 

The paper is organized as follows. In Sec. 2, we describe the models under consideration. In Sec. 3 we 
discuss the effects of the tree level operator Eq. (1.1) in neutrino physics. In Sec. 4 we discuss the effects 
of the operators induced at one-loop in charged lepton processes. In Sec. 5 we summarize our results. 

2 Right-handed neutrinos in extra dimensions 

We will study models with large fiat extra dimensions [1, 3, 4, 5] and models with one warped extra 
dimension [2, 6]. 



Large flat extra dimensions 

We consider (4 + (5)-dimensional massless fermions \I'j(x^,y) which, inside their 2^^"*" -''^ components (for 
even 6) or 2^^"'"")'^ components (for odd 6), contain the degrees of freedom of the "right-handed neutrinos" 
^ffi (^ = 1; 2, 3 is the generation index). The fermions ^i interact in our brane, through their components 
VRi, with the standard left-handed lepton doublet Li in a way that conserves total lepton number. The 
relevant part of the action is 



S= (fx d% *iiiZ)^i + 



d^x 



LiidLi 



LiXijURjH'^ + h.c. 



(2.1) 



where i3 is a 4 + 5 dimensional Dirac operator, H is the SM Higgs doublet in four dimensions and A 
is a matrix of Yukawa couplings with dimensions (mass)" . As is manifest from Eq. (2.1), A can be 
made diagonal without loss of generality at the price of introducing the usual unitary matrix [/, which 
describes flavour-changing charged current neutrino interactions. 
Using the KK decomposition 



^i{x,y) 
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yx) exp 



m ■ y 
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(2.2) 



where Vs = {2t:R) is the compactified volume, and performing the d y integration, Eq. (2.1) yields the 
four dimensional neutrino Lagrangian 



C = LiidLi + ^ I ^, 



n ■ 7 






A, 



V 



LiiyRHjH^ + li-c. 



(2.3) 



where 7 are the extra-dimensional Dirac matrices. After electroweak symmetry breaking, the neutrinos 
obtain a Dirac mass matrix rrii, = Xv/^/Vs, where v = 174 GeV is the vacuum expectation value of the 
Higgs boson. The heavy modes give corrections suppressed by Rm^, that are negligible in the cases of 
interest. 

In order to compute the tree-level effects due to the presence of the tower of KK states, we construct 
an effective Lagrangian by substituting in Eq. (2.3) the solutions of the equations of motion for the heavy 
fields, 

-1 A* 



^ r 



n ■ 7 



^LjH. 



(2.4) 



With an abuse of notation, we have identified Lj in Eq. (2.4) with a higher-dimensional spinor, in which 
Lj fills the components corresponding to the right-handed neutrinos, while all other components are zero. 
Ignoring higher derivative terms and summing over all KK states up to an ultraviolet cut-off A, we obtain 
the effective operator in Eq. (1.1) with coefficients 



(AAt),,A'5- 



2^2 



tjj 



Ss 

(2^ 



(2.5) 



where Ss = 27r '^/r(5/2) is the surface of a unit-radius 5-dimensional sphere, A is the dimensionful 
matrix of Yukawa couplings and A parameterizes some ultraviolet cut-off. We take A to be the KK mass 
at which we cut off the summation. The operator in Eq. (1.1) is gauge invariant, in spite of containing an 
ordinary derivative (rather than a covariant one), because HL is a SM gauge singlet. The coefficients eij 
of this operator have the same flavour structure as the physical neutrino masses myml oc AA^, although 
its overall factor is model-dependent. The order of magnitude of this single free parameter is fixed if — 
motivated by the hierarchy problem — we assume that all extra-dimensional physics is at the TeV scale, 
A~^' ^ K ^ TeV. In this case the coefficients eij in Eq. (2.5) are of order ls{v/ TeV)^, where Is is a typical 



loop factor in 5 dimensions. It is useful to write the coefficient of Eq. (2.5) as Cij = e{'m,^ml)ij / Am: 
because e is the only free parameter of the model (at tree level).* Explicitly, 
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(2.6) 



We now recall how this model can be linked with the Higgs mass hierarchy problem. 

If gravity propagates in d extra dimen- 
sions {d > 6, where all extra dimensions 
have the same radius R and the topology 
of a torus), the reduced Planck mass M4 = 
2.4 X 10^^ GeV is related to the reduced 
Planck mass Md in Z? = 4 + d dimensions 
by _ 

Mj = {2ttR)'^mI+'^. (2.7) 

It is useful to rewrite the dimensionful Yuka- 
wa coupling A in terms of a dimensionless 
Yukawa parameter A as 



A = A/M?/'. 



(21 



With this definition, the neutrino mass be- 
comes 



m. 



Xv 



M 



S/d 



D 



Ma 



(2.9) 



The simplest case 5 = d requires very large 
A or Md in order to obtain neutrino masses 
that satisfy the atmospheric neutrino data 
(see Fig. 1). On the other hand, for 6 = 5 
and d = 6, 



> 




1/2 . 


10^ 


r 


3/5.--^: 


10 2 


r ,.. 


--'"" _2/>^ 


10 


r .. 


■-'■'" 3^ 


1 


r 


^y^ 


-1 
10 


Y 


/^\:,'^\-xC^X^O<X^ 


-2 

10 


f^ 


1, 


-3 
10 


'■ 


/ : 


-4 

10 


- 


,.-''' - 


-5 
in 





■m,. 



A 



^)''3xlO-2eV, (2.10) 



10 



JO 10 

Md (TeV) 




10 10 
Mq (TeV) 



and the atmosph_eric mass scale is obtained Figure 1: m^/A (le^t) and e (right) as a function of Md for 
with A ~ 1 and Md ~ TeV, while the solar different values of S/d, assuming A = Md- The horizon- 
neutrino puzzle can be solved for values of ^Z band in the left panel shows the mass range selected by 
A slightly smaller than one. Figure 1 (left) atmospheric neutrino data /or A = 1. 

depicts values of rriy/X as a function of M/j, for different values of 5/d. Supernova bounds [8] force all KK 
neutrino states to be heavier than the typical supernova temperature T ~ 30MeV. In the present model, 
these bounds imply d ^ 3. The reduced Planck mass Mq is related to the phenomenological parameter 
M/5 used to study graviton effects at colliders as M^ = (27r) /'^"^^M/). The present collider bound on 
Md is Md^ 200 GeV (for d ~ 6) while the LHC can improve it by a factor ^10 [16]. Fig. 1 (right) 
depicts the values of e, defined in Eq. (2.6), as a function of Md in the case A = Md for different values 
of 5 and d, and for Aml^^ = 3 x 10~^ eV^. It is important to keep in mind that e is strongly enhanced 
if A > Md (it scales like (A/Md)^)- For a fixed ratio A/Md, e decouples like l/A^, as the new-physics 
scale increases. This behavior is clearly visible in Fig. 1. 

One may consider a few variations to the "minimal" model described above: 

(A) The right-handed neutrinos could have some extra dimensional mass term m > rui, [4] . The mass 
term m does not affect the dimension-six operators, but affects neutrino masses — now related to m 



'Throughout the paper, we will assume that the neutrino masses obey a normal hierarchy, i.e., mi <^ ni2 <^ ma, so that 
Am^tm ss mj and Aml^^ ^ ml. 



and A by a higher-dimensional "see-saw" relation, rrii, oc }? jvn. In these models, the e^j coefficients 
are not directly related to neutrino masses, and therefore contain additional mixing angles and CP 
violating phases beyond the ones in the neutrino mixing matrix J 

(B) Different massless 'right-handed neutrinos' could have a different UV cut-off, or live in a different 
number of extra dimensions. An interesting case is obtained with 2 right-handed neutrinos living, 
respectively, in 5 and 6 extra dimensions with equal radii. In this case one can reproduce the small- 
ness of the solar Am,^ with respect to the atmospheric Am,^ using comparable Yukawa couplings 
(see Fig. 1). This model contains 4 mixing angles and 2 CP-violating phases (rather than the 3 
mixing angles and 1 CP phase of the minimal model). The reason is that, unlike in the minimal 
model, it is not possible to perform flavour rotations of the right-handed neutrinos. 

Qualitatively, the main new feature of non minimal models is that the tree level exchange of bulk neutrinos 
again yields the operator Eq. (1.1), but the overall coefficient is now different for the "atmospheric" and 
"solar" contributions. Assuming that all Yukawa couplings are of order of the "fundamental" mass scale, 
in case (B) the "solar" contribution becomes comparable to the "atmospheric" contribution, i.e., all 
€ij in Eq. (1.1) are expected to be of the same order, contrary to the minimal model where e^T-je^p, oc 
^"^atm/^"^sun' assumiug that |C/e3| is negligible (this will be discussed in detail in the next section). 
Warped models generically give effects qualitatively similar to case (B), as we describe below. 

Warped extra dimensions 

A strong suppression of gravity could be generated by an extra-dimensional red-shift factor [17, 2]. We 
will concentrate on the simplest scenario [2], containing one extra dimension with topology S1/Z2 {i.e., 
a segment), parameterized by a coordinate y. The SM fields live on a brane fixed at one of its borders 
(y = R), and another brane lives at the opposite border {y = 0). After an appropriate fine-tuning of 
the cosmological constant A in the extra dimension and the tension of the two branes, the background 
metric has the form ds^ = e~ '^rj^ydx^dx'^ + dy"^ , where the warping constant k is determined by the 
five-dimensional cosmological constant A and gravitational scale M as /c^ = — A/24M^ [2]. By rescaling 
the metrics at 7/ = to its canonical form, one finds that all mass scales in the four-dimensional SM 
Lagrangian get red-shifted by a factor e~^^ (including the ones that should suppress unwanted non 
renormalizable operators), and that the four-dimensional Planck mass is Mpj = (1 — e~'^^^)M^/k. One 
can try to explain the gauge/gravitation hierarchy by stabilizing the size of the extra dimensions R [18] 
such that e-^^ ~ TeY/Mpi.^ 

It is again interesting to consider five-dimensional 'right-handed neutrinos' with Yukawa couplings A 
to the SM fields living on a brane. In this case it is necessary to give some five-dimensional Dirac mass 
terms m to the extra-dimensional neutrinos in order to naturally explain the smallness of the neutrino 
masses [6]. Despite this higher dimensional mass term, there is still one very light KK mode, and its 
effective Yukawa coupling with the active neutrinos is strongly suppressed if m > k/2. Then, the Dirac 
neutrino masses depend very strongly on m, rriy ~ ve~'^ , and neutrino masses which span many orders 
of magnitude are easily obtained even if all the higher dimensional Yukawa couplings are comparable. 

The other KK states have TeV-scale masses and unsuppressed four-dimensional Yukawa couplings to 
SM fields. They generate the operator Eq. (1.1) with coefficient 

e.i ~ ^-j—{\X%. (2.11) 

Because this is a 4+1 dimensional theory, the infinite sum over all KK modes is finite, and the introduction 
of an arbitrary ultraviolet cut-off is unnecessary. On the other hand, it is not useful to perform a precise 



^If m, '^ m„, more than three mass eigenstates participate in osciUations. Explaining the origin of a higher-dimensional 
mass term comparable to the neutrino masses is, perhaps, the biggest challenge for this type of scenario. 

'It has been conjectured that this model is equivalent to walking composite technicolour [19]. This dual version is not as 
aestetichally appealing, but its problems with experimental data are better known. 



computation because, similarly to the non minimal flat models, there is no direct connection between eij 
and the neutrino masses. The main conclusion is that all e^j have comparable values, e^j ~ 10~ , if A, 
k, M have comparable mass scales. Smaller Cij can be obtained for smaller Yukawa couplings A. 

We do not consider "non minimal" warped models, because they do not seem to lead to new interesting 
effects in neutrino physics.^ 

3 Tree level effects in neutrinos 

The tree-level operator Eq. (1.1) only affects neutrinos.* For this reason, Eq. (1.1) can give rise to 
detectable effects in neutrinos which are compatible with bounds from charged-lepton processes, affected 
only at the one loop order. After electroweak symmetry breaking, the operator in Eq. (1.1) contributes 
to the kinetic term of the neutrinos. The effective Lagrangian for the 3 left-handed neutrinos z/ and for 
the zero mode right-handed neutrinos Vfi becomes 

C = Vii&^Sij + eij)vj + VRiBuR - {vmyVR + h.c.) + ( —j=IVl/u + h.c. ) + -—vfv. (3.1) 

V V 2 / 2 cos Pw 

where i = {gl, IJ'LjTl) are the left-handed charged leptons in the mass eigenbasis, rrijy is the neutrino mass 
matrix, and eij are adimensional numbers defined in Eq. (2.5), expected to be of order ig{v/TeY)'^. The 
kinetic and mass terms can be simultaneously diagonalized by the following non- unitary field redefinitions 

i^R ^ Uri^r, ul ^ UlKU'^vl. (3.2) 

Here U are unitary matrices (such that Ul diagonalizes eij) and K is a diagonal matrix whose elements 
are 1/a/1 + e^, where ej are the eigenvalues of eij. In this new basis the neutrino kinetic and mass terms 
are fiavour diagonal, but the neutrino interactions with the gauge bosons are modified 

C = Didu + VRiduR - {pnmnyRn + h.c.) + ( —=lVwW^ + li-c- ) + t; j^^yz^^-, (3.3) 

V V 2 / 2 cos t^w 

where nin are the mass eigenvalues given by mvn/\/^ + ^n, Vw = UlKU'j^ and Vz = U'^K'^U']^. It is now 
easy to compute oscillation probabilities. The transition probability for a neutrino both produced and 
detected via a charged-current W interaction is given by 

2 



Pii^. ^ ^/) 



WTrW*JEnt 



Y.^^Zyfu^ 



f\rr) 
E^-Err. = ^^. (3.4) 



Differently from the standard case, this expression cannot be simplified to the usual form, because V^ 
is not a unitary matrix. 

The minimal model described in detail in the previous section predicts eij oc {mym\,)ij, such that 
C/^ = E. The transition probability reduces to the more transparent form 



P{n ^ vf) 



TT- TT* 
y^ ^mUfn ^ ^iE„t 



+ e?7i|/Am^ 



atm 



(3.5) 



nn more general extra dimensional black-hole backgrounds, the space warping factor can be different from the time 
warping factor [20]. This could lead to right-handed neutrinos that travel with a velocity c(l -I- 5c) different from light. 
This model has a conjectured holographic dual [19] where Lorentz invariance is broken in a more obvious and old way [21]: 
normal right-handed neutrinos have an index of refraction n due to interactions with an 'ffither' (composed of some hot 
conformal matter). At tree level and up to negligible m/E corrections, these effects do not affect neutrino oscillations as 
long as there are only Dirac mass terms [22], or more generically as long as the energy eigenstates do not mix left with 
right-handed neutrinos. Even in models where this not the case, neutrino effects caused by a non universality of the speed 
of the light do not seem phenomenologically interesting because have an energy dependence (different from the observed 
one [23]) that allows to derive very strong constrains. Upward through going atmospheric muons in SuperKamiokande 
(with L r^ 10'* km and Ei, ~ TeV) put the bound 5c%,l/LEi, ~ io~"*~^^' (for large mixing between flavour and velocity 
eigenstates), preventing detectable effects in planned new neutrino experiments. 

*The operator in Eq. (1.1) does not contribute to the h ^ ui? decay, due to conservation of angular momentum. Only 
the emission of KK states lighter than mn yields h -^ vun corrections to Higgs decays [5, 24]. 



where U is the usual unitary neutrino mixing matrix. Notice that the summed survival probability 
J2f=e,iJ.,T ^i'^i "*■ ^/) — 1 — 2e is not equal to 1, because of the non-unitary e effects. This corresponds 
to a non- vanishing transition probability into sterile KK modes or, in the effective theory, to a reduced 
interaction of the neutrinos with the W boson. 

Next, we discuss a few possible experimental signals, concentrating on searches for neutrino oscilla- 
tions. 

Flavour transitions at very short baselines 

In terms of the Cij parameters in Eq. (1.1), the transition probability at a very short baseline (L ~ 0) 
derived from Eq. (3.4), for small €ij, is given by 

P{i^i -^ Uj) ~ |eij|^ (i 7^ j), P{vi ^ i/j) ~ 1 - 2eii (3.6) 

The present bounds on the flavour violating eij = e^j, coming from neutrino experiments, are 

|e^^| < 0.013, \eer\ < 0.09, \ee^\ < 0.05. (3.7) 

The dominant bounds on e^r and egr are due to the Nomad experiment ^^ [25], P{v^ — > i't) < 1-6 x 10~^ 
and P(fe -^ Vt) < 7.4 x 10~^ at 90% CL. The bound on ee^ comes from Karmen [26] and Nomad [25]. 
The LSND anomaly [13], which may be interpreted as evidence for P{h'^ -^ i/g) ~ 2.5 x 10^'^ could 
tentatively be solved if |ee^| ~ 0.05. This is in slight conflict with the bound from Karmen, and in 
strong contradiction with the bounds from jjl — > 67, which we will discuss in the next section. 

The diagonal elements en = e*j are constrained by Chooz [27] and Bugey [28], yielding ege < 0.025 
and, because of the modified neutrino couplings to the Z boson, by the invisible Z width measured at 
LEP [29], which leads to \eee + e^^ + e^^| < 0.013 at 90%CLt. The liWv couplings are a factor 1 - eii/2 
smaller than in the SM, so that lepton universality tests [30] in r and vr decays give 



^^l^l 



0.0027 ± 0.0024, e^^ - e^r = 0.0003 ± 0.0042, eee - e^r = 0.0011 ± 0.0045. (3.8) 



A slightly more stringent bound, egg + e^^ < 0.002, can be obtained by comparing a global fit to the 
precision LEP data with the /i lifetime. The bounds on the flavour violating eij from Z, r, /x, vr decays are 
not competitive with the ones from neutrino oscillation searches. 

In the flat minimal model, the e^j parameters are determined in terms of the single unknown parameter 
e, defined in Eq. (2.6), and of the measurable (and already "partially" measured) neutrino oscillation 
parameters as 







[C/gsl e-^'^Am2^/ 



(3.9) 
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where we assumed maximal solar and atmospheric mixing. Here C/gs is the element of the neutrino mixing 
matrix which is currently constrained to be small by the Chooz reactor neutrino data [27] and (p is the 
CP-violating phase. The bounds previously discussed can be turned into constraints on e, in the case of 
the minimal model. It turns out that the most stringent bounds are e < 0.004 (from /i, r, vr decays and 
precision LEP data) and e^0.03 (from neutrino experiments). 



'The CERN neutrino experiments, Nomad and CllORUS, were motivated by theoretical prejudices for small mixing angles 
and for warm dark matter. Today they are among the most significant probes of extra-dimensional neutrinos. 

"* Since the present measurement of neutrino counting at LEP agrees with the SM only at the 2-cr level, and a nonzero 
e reduces the effective number of neutrinos to Ni, = 3 — 2 Tr tij , this measurement gives a possible indication of a positive 
effect, fee + e„„ + e^r = 0.008 ± 0.004. 



The most sensitive future experimental search seems to be r or wrong signed /i appearance at a future 
near-detector of a i/-factory. Such a detector, located at < 1 km from the neutrino source, has already 
been studied in detail as a tool for neutrino deep inelastic scattering experiments [31]. It expects to 
observe 10^ charged current SM events per year. It seems possible [31] to probe e^r and e^r values down 
to the level of 10~^, by looking for r appearance, and ee^ down to similar values by looking for wrong 
sign muons, therefore improving the present bounds by two orders of magnitude. More dedicated studies, 
however, are still required. The sensitivity on the flavour diagonal e^^ from searches in the disappearence 
channel seem to be slightly worse, being limited by the uncertainty on the incoming neutrino flux. The 
sensitivity on e^e should be weaker than the one to e^^, while e^r effects cannot be probed. 

CP violating effects 

In order to understand qualitatively what the general exact formula Eq. (3.4) means in the various 
possible models, it is useful to specialize it to the short-baseline case 5ij <^ 1 where 5ij = AmfjL/2Ey. 
Furthermore, it is useful to count how many CP-violating phases are contained in the €ij = e^j parameters. 

1. In processes where neutrino masses can be neglected, the €ij contain {Ng — l){Ng — 2)/2 CP-violating 
phases (where Ng = 3 is the number of generations) which do not give rise to any observable effects. 

2. If the eij can be neglected (normal oscillations), the neutrino mixing matrix contains the usual 
{Ng — l){Ng — 2)/2 CP phases. For Ng = 3 there is one CP-violating phase, (f), and the CP- 
asymmetric part of the oscillation probability is the same in all flavour transitions fj — > Vj ^ fj, 
and, for small L, is proportional to L^: 

_ P{i>i -^ Uj) - P{vi -^ Uj) . Si2 . 6i3 . 623 , . 

PqP = '—— '— = 8Jcp sm — sm — sm — « JcpOi20i3()23, (3.10) 

where Jqp ~ —-^UessiiKp (assuming maximal mixing in atmospheric and solar oscillations, and 

\Ue3\ « 1). 

3. In a generic process affected by neutrino masses and by the e^j, the neutrino mixing matrix contains 
the usual {Ng — l){Ng — 2)/2 CP phases, and the eij parameters contain other Ng{Ng — l)/2 CP 
phases. Therefore, unlike the standard case, CP violation can be present when only the dominant 
two-generation z/^ <^ v^- atmospheric oscillations is "turned-on." 

For example, ignoring subleading (5sun effects and keeping terms up to second order in (5atm 



P{^, 



sm26'atm^ 



, ^ = 10"^ A..^ LgOGeV^ 

2 3 • 10-3 eV^ km ^^ ^ ' 

The fact that these two effects could be comparable offers an opportunity to measure CP-violating effects 
present, if e^,- is complex, as a difference in P{i'^ — > ^t) from P{i'f_i -^ i>^) proportional to L. The CP 
asymmetry can be maximal with a suitable choice of the pathlength and neutrino energy, if e^r has a 
large CP-violating phase [32]. 

In the minimal model presented in the previous section, the only CP-odd phase is the one present in 
the neutrino mass matrix (recall that e^j oc {m,yml)ij), and new CP-violating effects (proportional to L) 
are suppressed by Jqp and by two powers of e: 

PS3P = Jcp[-25i2e2 + (5i25i3<523]. (3.12) 

Given the current constraint on e<0.01, such effects are negligible. 

For specific non minimal models, it is necessary to check whether the (indirect) relation of Cij and 
the neutrino mass matrix allows to "rotate away" potentially physical phases in eij. This is the case if 
a single right-handed neutrino gives a dominant contribution to e^, as happens in the minimal model. 
Furthermore, in order to obtain observable effects, it is necessary to generate a large eij with a large 
phase, while keeping |ee^| small, due to the severe constraints already imposed by searches for /i — > 67 
(see next section). In spite of all the difficulties, however, it is possible to build specific models which 
yield significant CP-violating effects. 



Matter effects 

Non-flavour diagonal interactions with the Z, W bosons give rise to non standard MSW corrections. 
The corresponding "matter potential," however, is too small to solve the solar neutrino puzzle or affect 
significantly the SM effect due to charged current i^e — e scattering. On the other hand, the SM matter 
effects in i/^ <-^ v^- oscillations are suppressed by ~ [m-r/TrMiY)'^ [33], so that non SM effects could be 
dominant and perhaps detectable. In the minimal flat model, f^ ^^ Ur oscillations in neutral matter are 
modified to 



P(i/^ -^ Uf,) = 1- sin (26'atm) sin 



'-S^--,/2GfNAl 



AE^ o ■ - . - e 



+ 0{e'), (3.13) 



and could result in a narrow resonance in long baseline oscillations of high energy anti-neutrinos, E^, ~ 
30GeV/e. Similar effects have been discussed in [34]. 

Matter effects are also known to substantially change the expected neutrino fluxes from super- 
novae [35] . The new matter effects due to the KK right-handed neutrinos produce a new MSW resonance 
in supernovse which seems, however, not to perturb the observable V(, and u^. spectra in a significant way. 

4 One-loop effects in charged lepton processes 

The effects of any high energy theory can be described at much lower energies by higher dimensional 
operators. The relevant set of operators and their coefficients depend on the ultraviolet cut-off which is 
used to make computations in the low energy effective theory. Using dimensional regularization (so that 
no power divergences arise in loop computations) , the low energy effects of bulk neutrinos are described 
by a few dimension-six operators, including iilj^f magnetic moment operators, iiijZ vertices, and four 
fermion interactions. 

These effects are partly due to computable SM loop corrections to the operator in Eq. (1.1), gener- 
ated at tree level. However, comparable contributions come from high energy effects in the full theory 
containing the KK modes of the higher dimensional fermions. As before, it is important to emphasize 
that the coefficients of the effective operators cannot be computed, since we do not have a renormalizable 
high-energy theory. The coefEcient of each operator is a free parameter that cannot be calculated in 
terms of the non-renormalizable Yukawa couplings A of the extra-dimensional neutrino. At best, they 
can be estimated by introducing some explicit cut-off A. It proves particularly useful to rewrite the 
dimensionful Yukawa coupling A as 

because we can interpret the dimensionless couplings A as parameters for understanding how "strongly 
coupled" the right-handed neutrinos are.* As defined in Eq. (2.5), is is a (5-dimensional loop factor (for 
example £4 ~ l/(47r)^). All virtual effects are comparable to the effects of a four-dimensional heavy 
right-handed neutrino with a Yukawa coupling A (so that A ~ 47r corresponds to strong coupling). The 
effective Lagrangian at energies smaller than the cut-off A is 

C ~ ^{LH*dLH) + ^ {eX^ifi ^ 67] + ^^(A^ + g^) [gv^[Z^xe] + [fi ^ 3e]) + ~X^g^[fiq ^ eq]} . (4.2) 

We have used a shorthand notation to indicate the various SU(2)i (8) U(l)y-invariant operators. For 
example, gv'^[Zfie] indicates operators like {H*D^H){L^^L). While the magnetic operator [fi — > 67] is 
generated at order A^, other dimension-six operators also receive contributions at order A^. In the case 
of the minimal model presented in Sec. 2, the loop effects are explicitly calculated (making use of a hard 
ultraviolet cut-off) in the Appendix. 



*When in section 2 we discussed the connection of neutrino masses with gravity, it was convenient to parameterize the 
dimensionful fundamental Yukawa couphng A in terms of another dimensionless coupling A, defined using the reduced Planck 
mass as the unit of mass. 



Note that all dimension-six operators decouple as 1/A^ when A -^ co. The correct decoupling is also 
found in the expressions presented in ref. [15] , once the definition of the mixing parameters is taken into 
account. 

We remark that we are only considering the effects due to the extra-dimensional right-handed neutri- 
nos: a full quantum gravity theory at TeV energies is expected to give additional effects. In particular, 
it is important to note that while some unknown mechanism could be responsible for suppressing effects 
which lead, e.g., to proton decay, it is hard to believe that there are no extra contributions to lepton 
flavour violating processes, given that the neutrino data indicates that individual lepton flavours are 
strongly broken symmetries (see [36] for a discussion of these effects and possible ways of addressing this 
issue). Nonetheless, we will ignore here such contributions, which are impossible to estimate. 

In spite of all the intrinsic uncertainties, useful results can be obtained. Both the tree level and 
the more relevant one-loop operators are suppressed by the same order of A/A. This implies that one- 
loop effects are only suppressed by a factor of order l/lGvr^ with respect to the tree level term.T Some 
operators, like the four-fermion operators contributing to ;U — > 3e and /i — > e conversion in nuclei, may in 
fact be enhanced with respect to ;U — > 67 by X'^/g^ if A > (7. Moreover, the coefficients of the operators 
included in [Z^e] and [/j, —>■ 3e] receive a lnA'^/m^^ enhancement. Extra-dimensional models predict (up 
to order one factors) relations between flavour violating effects in the charged and neutral lepton sectors. 
In the simplest "flat" model we discussed in Sec. 2, predictions for many charged lepton flavour violating 
processes will be strictly related to the observed neutrino oscillation parameters (Am^, and neutrino 
mixing angles), such that, e.g., it is possible to predict BR(r — > fj,'y)/BK{fi — > ej). 

Therefore, minimal extra-dimensional models are more predictive than other beyond-the-SM sources 
of lepton flavour violating phenomena related to right handed neutrinos. In the constrained MSSM, 
for example, the presence of very heavy right-handed neutrinos (that generate small neutrino masses 
via the see-saw mechanism) yields potentially large flavour violating effects in the muon sector. While 
all branching ratios are precisely calculable in terms of the parameters of the theory (the MSSM is a 
perturbative gauge field theory), it is impossible to establish a connection between the observed flavour- 
violating neutrino masses with predictions for fj, — > ej, etc, due to the presence of too many unknown 
flavour-mixing parameters in the right-handed neutrino sector. 

In this section, we concentrate on charged lepton flavour violating phenomena, and also comment on 
the anomalous magnetic moment of the muon. Note, however, that the effective Lagrangian Eq. (4.2) also 
allows e~^e~ — > if^J {i,j = e, /i,r) processes [15]. Such effects could be studied at the Z resonance with 
a next-generation e^e~ collider. However, searches for processes with i ^ j do not seem very promising 
given the current bounds on Cij [37], and the effects on processes with i = j seem less sensitive than 
corrections to Z observables, which are generated at tree level. 

Explicit expressions for the processes of interest are derived and listed in the Appendix. They agree, 
apart from minor differences, with the corresponding expressions in [15] when the same model and 
the same arbitrary cut-off is chosen. Our numerical results are, instead, different. In particular, the 
experimental bounds we obtain are much weaker and do not, in general, require fundamental scales 
above 100 TeV or 10 TeV. We discuss some of our results in what follows. 

Flat extra dimensions 

We will concentrate on the minimal model outlined in Sec. 2. This model contains only two new free 
parameters, e and A. Up to order one factors, the rates for the different charged lepton process depend 
only on these two parameters, on the neutrino masses, and on the elements of the standard neutrino 
mixing matrix Uai. The largest Yukawa coupling is then determined by the relation 



-X = ^^e{6-2). (4.3) 



'This nice property is not shared by virtual effects mediated by higher dimensional gravitons, that generate dimension 6 
operators at loop level but only dimension 8 operators at tree level. 
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As mentioned before, we will assume that the neutrino masses are hierarchical (mf <C m^ <C m^) such 
that all information regarding neutrino parameters can be obtained from neutrino oscillation experiments. 
We later comment on non minimal flat models and warped models. 

/x — > 67 and r — > //7 

Using the expressions derived in the Appendix, the branching ratio for ^u — > 67 can be written as follows: 



BR(/i-e7) = pe' 

Svr 



EUejU^ 



m^ 



M ^ 



m 



atm 



3a 

8^* 



^ Am 






atm 



(4.4) 



Similarly, the branching ratio for r -^ /_f7 is given by 



BR(r -^ fi-y) 



0.174^.2 

Svr 



T.u,,u, 



m^ 



TJ 



\m. 



atm 



0.174^6^ 

Svr 



U,2U:^ 



^"^smi I jT jT* 



\m. 



atm 



(4.5) 



The branching ratio for /_f -^ 67 depends on the solution to the solar neutrino puzzle (i.e. Ue2 and Am^^^) 
and on the unknown C/gS element. On the other hand, BR(r — > /i7) depends dominantly on the well 
measured atmospheric parameters:^ l^/xst^-sl ~ 0.5 and Am^^^^ ~ 3 x 10~^eV^, such that 



BR(r -^ fi-f) ~ 0.174— e^ PusK^f = 4 x 10"^ e^. 
Svr 



(4.6) 



The bound Eq. (3.7) from neutrino experiments on e^r ~ e/2 implies, in a model-independent way, that 
a possible effect in r — > /X7 is at least two orders of magnitude below the present limit BR(t -^ fi'y) < 
1.1 X 10"^. 

The dependence on the free parameter e cancels out in the ratio between BR(/_f — > ej) and BR(r — > 
/X7), that is expressed only in terms of neutrino oscillation parameters. 



BR(/i -^ 67) 
BR(r -^ ^7) 



Ue2U;^Aml^ + Ue3U;,Aml,^ 



0.17A\U.3Ur3Aml 



(4.7) 



atm I 



The numerical value depends on the still unmeasured Ues and on which is the true solution of the 
solar neutrino puzzle. The dominant contribution to /x ^ 67 could be related to the "atmospheric" or 
"solar" neutrino mass splitting. If the solar anomaly is due to LMA oscillations, and if |C/e3| is such that 
Amg^jj ^ [/es-^^atm' the "solar" contribution dominates, and we predict 



BR(/x ^e-f)^lx 10~^e^ 



BRifi ^ ej) ^ ^ ^ ^^_^ ( \Ue2U^2\ 0.5 



^»^sun 



3 X 10-5 eV 

, 2 / 



^"^sun 



3 X 10-3 eV^ 



^"^atm 



(4.^ 



BR(r^/x7) \^ 0.35 {U^^Ursl J V^ x lO^^ eV^ 

When compared to the current experimental bound BR(;U — > 67) < 1.2 x 10"^^, we obtain the limit 



(4.9) 



e < 3 X 10~2 



'3 X 10-5 eV^ 

A?n2^„ 



II \Ue3\ < -T 2 • 



(4.10) 



In the near future the experimental sensitivity to BR(;U — > ej) is expected to reach 10 [3S]. 

^ The r — > 67 branching ratio is suppressed by the unknown reactor angle Ues, so that it is smaUer and more uncertain 
than the r — > /i7 branching ratio. 
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On the other hand, if the solar parameters fall in the LOW or SMA regions,'^ and \Ue3\ is large enough, 
the "atmospheric" contribution to BR(/i — > 67) dominates, so that 



BR(,^e7)«4xlO-e^|^e3/0.3P, |||^^ 



Ues/O.S 2 



Ur3/0.7 



(4.11) 



In this case the bound on e becomes 

e<5xlO-^|0.3/C/e3| if |C/e3| » AmL/Am^,^. (4.12) 

If Ue3 is close to its current experimental upper bound, \Ue3\ ^0.3, the bound on e is so strong that 
experimental signals in the neutrino sector can only be very small. 

Fig. 2 shows the values of BR(/i -^ 67) and BR(r — > fi'y) as a function of {Uesl ^^^ Amgy^/Am^^„^ = 
10~^, iL^/is/C/rSp = I5 \Ue2/Uei\'^ = 2/3 (as Suggested by atmospheric and LMA solar oscillations), and 
no CP violation in the neutrino mixing matrix, for different values of e (and A). If the "solar" and 
"atmospheric" contributions to // ^ 67 are comparable (which happens at Ues ~ 10~^^~^' for LMA solar 
oscillations), the CP-violating phase in the neutrino mixing matrix may either enhance or suppress the 
branching ratio, depending on how the two terms interfere. 

// — > eee and fi ^ e conversion in nuclei 

Unlike the ii — > ij-y decays, these rare leptonic processes are more dependent on the unknown ultraviolet 
details of the models, but still a few interesting results can be obtained. 

In // ^ eee process, several operators may contribute significantly (see the Appendix for detailed 
expressions). First of all, the // — > e7 magnetic penguin operator contribution to /i — > eee is enhanced 
by ln(m^/me), which is a consequence of a collinear divergence of the electron-positron pair in the 
me — > limit. Second, the Z and photon penguin diagrams are enhanced by an ultraviolet divergence 
~ ln(A/r?H4/). When A is significantly larger than mw, this log-enhancement is important. Finally, as 
discussed after Eq. (4.2), the Z-penguin and box contributions contain A^-terms. If the higher-dimensional 
KK neutrinos are coupled strongly enough such that X > g, then the A^-terms may dominate over the 
other contributions to BR(/i — > eee). 

In light of this discussion, the ratio of branching ratios BR(/i — > eee)/BR(;U — > e7), in which the 
overall e^ dependence in both processes cancels out, can be written as 

BR(/x -^ eee) _ ^ , ^ , ^ / . , on 

Here, C^^e7 is the contribution from the /x — > e7 magnetic penguin operator 

The second term CinA contains the InA/niw contributions from the Z and photon penguin diagrams of 
// — > eee. For example, the Z-penguin diagram leads to 

ClnA ^-3— 3- _ ^ ' n ■ 2 

8tt \ sm^ 9w 2 sm^ 

We find that CinA = 0.03 (0.2) for 6 = 5 and A = 1 TeV (10 TeV), when all In terms present in the Z 
and photon penguin diagrams are included. 




^ After SNO [12], the SMA solution is strongly disfavoured by the solar neutrino data [39]. 
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Figure 2: The most interesting leptonic observables as a function of \Ue3\, for 5 = 5, e = 0.003 (top) 
or 0.0003 (bottom), and A = 1 TeV (left) or 10 TeV (right). We assume Am^^JAmlt^ = 10"^, 
lU^s/Ursl'^ = 1, \Ue2/Uei\'^ = 2/3 (i.e., maximal mixing in the atmospheric sector and the LMA solution 
to the solar neutrino puzzle), no CP-violation in the neutrino mixing matrix, and hierarchical neutrino 
masses (mf <C ttt-I ^ "^ij- Note that A is the largest Yukawa coupling, P^ut — > Vj) is the neutrino 
conversion probability at very short baselines, and the /i — > e conversion rate is computed for ^'^Al. 
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The last term in Eq.(4.13), C\, contains the contributions from A'^ terms in the Z-penguin and box 
diagrams: 



C.oc 



xV 



Ue2U;,{AmlJAml,J +U,iU;^ 



2- 



(4.16) 



When \Uez\ is small, C\ is small (suppressed by the small ratio of neutrino mass-squared differences 
squared). On the other hand, if \Uei\ is sufhciently large, C^ may in fact be the dominant contribution 
to /x ^ eee due to the potentially large iX/g)'^ enhancement. Numerically, C\ = 0.6 for A = 10 TeV, 
|C/e3| = 0.1, and e = 0.003, which corresponds to A = 5.4. 

Because of the contributions CinA and C\, the ratio BR(/i -^ eee)/BR(^ -^ 67) can be significantly 
larger than 10~^ if A ^ 1 TeV or A > 5. This is very different from predictions for lepton flavour violating 
processes from SUSY models with slepton flavour mixing, in which the contribution C^^e'y is almost 
always dominant [40] and therefore BR(^ — > eee)/BR(/i -^ 67) ~ 6 x 10"'^. Large contributions to 
BR(^ -^ eee) may be obtained in SUSY models with iZ-parity violation [41]. 

The fi ^ e conversion rate in nuclei behaves similarly to the fj, — > eee branching ratio: for A >> mw, 
the log-enhancement in Z and photon penguin diagrams is important and perhaps dominant, while if 
X > g, the A^ term in the Z-penguin and box contributions can be significant. In both fi -^ eee and 
/x — > e conversion in nuclei, the Z-penguin contribution tends to dominate because of the In A and A 
terms. When this is the case, the ratio R(/i — > e)/BR(/x — > eee) does not depend on the unknown 
ultraviolet details of the models. We verified numerically that the ratio is almost constant, varying in the 
range 10-13 in the case of /U — > e conversion in ^^Al, and 20-25 in ^^Ti, in a large region of the parameter 
space. This feature is a definite prediction of the models under consideration. Again, the situation is 
different from SUSY models with slepton flavor mixing, where one expects R(// -^ e)/BR(/i — > eee) < 1. 

Figure 2 shows the branching ratio of /i — > eee and the rate for ^ — > e conversion in ^^Al as a function of 
[C/esj for different values of A and e (we also fix (5 = 5). All of the features we have discussed can be readily 
observed. First, one can clearly see that R(/i — > e in ^^Al)/BR(/i — > eee) is roughly constant (~ 10) for 
all the depicted values of e and A. Second, at small values of [f/esj (where the A^ terms are negligible), 
the effect of the ln{A/mw) enhancement is visible: at A = 1 TeV, BR(^ -^ e^) >R(^ — > e in ^'^Al), 
while at A = 10 TeV the situation is reversed. Finally, the A^ enhancement is also clear if one compares 
R(/i — > e in ^'^A1)/BR(^ — > e7) and BR(/i -^ eee)/BR(/i — > e7) at small and large values of [C/esj, for A 
large. This behavior is more easily observed for e = 0.003 and A = 10 TeV. Indeed, we have verified that 
at even larger values of A, BR(/i — > eee) exceeds BR(/i — > e7). Note also that, for all values of e and 
A depicted in Fig. 2, the rate for // — > e conversion in •^^Al is larger than the proposed sensitivity reach 
of the MECO experiment [42], even in the limit of small [^esl (assuming the LMA solution to the solar 
neutrino puzzle). 

Finally, we point out that, in general, the "atmospheric" and "solar" contributions to /i — > eee 
have different CP-violating phases. Furthermore, their relative weight in the magnetic penguin operator 
is different from their relative weight in the four-fermion operators in Eq. (4.2), unless X <^ g. If 
the "atmospheric" and "solar" contributions are comparable, the interference between them produces 
observable CP-violating effects in polarized fj, — > eee decays [43] . 

g — 2 of the muon 

As derived in the Appendix (see also [44]), the contribution of the KK tower of neutrinos to the muon 
anomalous magnetic moment is 



Sa - ^^ '"^ rVlL^ P "^"^ - IQ-^c f4 17l 



''W 1 ^'"-atm 



Here we have assumed maximal mixing in the solar sector and neglected C/gs and Aml^^/ Aml^^^ correc- 
tions. Note that the sign of the effect is negative, in contrast to the BNL experimental result [45], that 

14 





present bounds from experiments with 


future sensitivity from experiments with 




neutrinos charged leptons 


neutrinos charged leptons 


kerl 
Kr\ 

€rr 


< 0.05 (Karmen) <, IQ'"^ {n -^ 67) 
<0.09 (NOMAD) <10-i (r ^ 67) 

< 0.013 (NOMAD) <10-i {t ^ fij) 

< 0.025 (reactors) < 10"^ (^ (j^ta) 

— ^ 10-3 {Z data) 

— <10-2 (Z^uu) 


~ 10-"^ {u factory) ~ lO"''^'' (/i decays) 
~ 10-^ (z^ factory) - lO^^-^ (^ ^ g^) 
~ 10-"^ (zv factory) ~ lO^^-^^ (^ ^ ^^) 
~ 10-3 [u factory) ~ 10-3-^ (^ data) 
~ 10-3 (zy factory) ~ 10-3-^ [z data) 
— ~ 10-3 (Z ^ Z/P) 



Table 1: Bounds and future sensitivity reaches for the eij coefficients defined by Eq. (1.1 



claims a discrepancy 6a^ = +(4.3 it 1.6) x 10"^ with respect to SM predictions (afflicted, however, by 
significant hadronic uncertainties). On the other hand, the present bounds on the Cij parameters already 
require in a model-independent way that the right-handed neutrino correction to (7 — 2 is smaller than 
the theoretical uncertainty in the hadronic SM contributions. If, in the future, more precise experimental 
and theoretical results establish the presence of a non SM correction to 5 — 2, extra dimensional mod- 
els could still account for it by invoking ad-hoc dimension-six operators like -r7(Lr"cr^i/-^ L)14^" with 
A ~ (3 + 4) TeV. There is no direct contradiction between such non renormalizable operators and bounds 
from the more sensitive precision data. 

Non minimal flat models and warped models 

Nonminimal flat models and models with warped extra dimensions do not allow one to relate the coeffi- 
cients of the dimension-six effective operators to the parameters in the neutrino Dirac mass matrix. For 
this reason, it is not possible to make interesting predictions. In non minimal flat models (B) and in 
warped models the naive expectation is that all eij are comparable, but is easy to avoid this conclusion. 
It is useful to estimate the ii — > ij^ decays in terms of the €ij parameters defined in Eq. (1.1). Up to 
order one uncomputable factors, ii - 



j7 experiments give the following bounds 



'^efi\ rO 



<io- 



^flT 



<io- 



(4.18) 



Such bounds prohibit the interpretation of the LSND anomaly [13] as due to Ce^, as discussed in the 
previous section. In fact, the present bound on Cefi is so strong that it will be hard to observe its effects 
in Ue ^-> I'fi transitions, even with a neutrino factory. 

5 Conclusions 

We have shown that the various minimal and non minimal models with right-handed neutrinos in flat or 
warped extra dimension are described at low energy by neutrino masses plus a specific set of dimension 
six operators. Up to order one factors (that are anyhow uncomputable since these models are not 
renormalizable), flavour conserving effects are described by three positive numbers €ee, ^fifi, ^tt, and 
flavour-violating effects are described by three complex numbers Ce^, Cer, e^r, deflned in Eq. (1.1). This 
gives rise to several relations between non-SM effects in neutrino observables and in lepton flavour- 
violating processes. Effects in charged leptons are suppressed by a one loop factor with respect to effects 
in neutrinos: 



P{i^^ 



Uj;L 



0) 



BR{£i -^ £j-f) ~ \eeij/47r\ 



{i^j) 



(5.1) 



Present bounds from neutrino experiments and from charged lepton processes are summarized in Table 1. 
Due to the loop factor, detectable neutrino effects are compatible with lepton flavour violating bounds, 
unlike what is obtained with a generic larger set of SU(2)j^-invariant dimension six operators, where 
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neutrino and charged lepton effects botli arise at tree level. Values of eij > 10~ (including possible 
CP-violating phases) will be probed by future neutrino experiments. The importance of ^ — > eee and 
H ^ e conversion relative to /i — > 67 can be enhanced, with respect to the usual magnetic-penguin 
dominance approximation, if the right-handed neutrino is strongly coupled or if the cut-off of the theory 
is significantly larger than the W mass. 

All these effects can be generated in four dimensions by adding "right-handed neutrinos" with TeV- 
scale masses and order one Yukawa couplings, but such choice of parameters is not motivated by neutrino 
masses. On the contrary, extra-dimensional models that try to address the hierarchy problem and 
to generate neutrino masses give an order-of-magnitude expectation for the eij parameters of eij ~ 
is{v/TeV)'^ ~ lO"^*^"" where £s is a loop factor in 6 dimensions. Therefore, these models are compatible 
with a "natural" value for the fundamental scale. 

In the flat minimal model, the six eij parameters are related and can be expressed in terms of a single 
unknown e, see Eq. (3.9). The value of e is at present constrained by LEP and neutrino experiments 
(e < 0.004) and by the /i ^ 67 decay, see Eqs. (4.10) and (4.12). Improvements in the sensitivity on rare 
muon processes and measurements at a future neutrino factory will significantly extend the probe on the 
hypothesis of an extra-dimensional origin of neutrino masses. 

Acknowledgments We thank Paolo Criminelli, Riccardo Rattazzi, Serguey Petcov and Alexei Smirnov 
for useful discussions. 

Note added The NuTeV collaboration [48] has recently reported a 3-cr anomaly in f^-nucleon scatter- 
ing, that can interpreted as due to a ratio between the v^Zv^ and the P^Wfi couplings (0.58 it 0.21)% 
lower than in the SM. Since LEP found that the Z couplings to other fermions agree with SM predictions 
at the per-mille level, the NuTeV anomaly could be due to physics beyond the SM that mainly affects 
the neutrinos. This is what happens in the extra dimensional models considered here, where the v^Zv^ 
coupling is 1 — e^^ smaller than in the SM, while the v^W fi coupling is smaller by only 1 — e^^/2 (see 
section 3). The NuTeV anomaly could be fitted by e^^ = 0.0116 it 0.0042. This value gives a reduction 
of the Z ^> vv width compatible with the LEP measurement, 2-0" lower than its SM prediction (see foot- 
note in page 7). However, the fact that the charged current is also modified imposes severe constraints 
(the strongest bound, e^^ -|- e^e % 0.002, comes from a comparision of the muon lifetime with precision 
electroweak data), which prevent a clean explaination the NuTeV anomaly. 

A Integrating out right-handed neutrinos 

In this Appendix, we present in detail the expressions for the effective operators that mediate rare muon 
processes and the muon anomalous magnetic moment when the entire tower of KK right-handed neutrinos 
(up to some arbitrary cut-off) is integrated out. The computation is done in the minimal flat model, but 
can be easily generalized to other cases of interest. 

Magnetic moment type operator for /i ^> 67 

The rare decay ^ — > 67 takes place at the one-loop level. The following magnetic moment operator is 
generated when all one-loop diagrams involving KK neutrinos are added. 

£ = ra^,eG^p{ALPL + ^ijPR)/iF"^ + h.c. (A.l) 



where Ai^ = 0, 



Ar = ^.^HyXG.imr/ml.), (A.2) 



647r2 ml, ^ ^" ^" 



^ , , 3;(1 — 6x -|- 3x^ -|- 2x^ — 6x^ Inx) ,. , 

G.(-) = -^ ^^yT^, '- (A.3) 
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Here the summation of all KK modes, which are labeled by the generation j and the integer vector n in 
the 5-extra dimension, are taken into account. m]j^ is the mass of a KK neutrino labeled by j and n, 
and Vji is defined as V/^ = UijVl"', where Uij is the standard active neutrino mixing matrix and Vl"' is 



,(i) 



the active-KK mixing matrix. When m^ ^ Tn,yj, m^ 
ratio for /x — > 67 is given by 



(J)2 



n' 



''/B? and V^" ~ m^j/m^ . The branching 



BR(/i -^ 67) 



487r2|Ai?|2 3a 






27r 



E^^G,(mg)Vm 



ly. 



nj 



(A.4) 



^r, V;^ - K';, Vl^ - ^i, and 



rn ' en 



The same expressions can be used for t ^> l'^ after replacing m^ 
BR(^ ^ euei^fj) -^ BR(t — > IuiVt)^ for Z = e or ^u. 

Next, we sum over all KK modes. In models with large extra dimensions, the sum over the KK states 
can be accurately replaced by an integral 



n|<AR 

E / 



n 



Jo 



(A.5) 



where / is any function and Ss = 27r '^/r((5/2). The cut off A is expected to be of order of the fundamental 
scale A ~ Md- Therefore, 



n,j j 

- T.UejU;^ml^SsR^ J dEE 



m 



i^3 



E 

<5-3 



G^{E?lml^) 



Ss 



25/7 



6-2^1110) \Md) A2(27r) 



Y.UejU;^mlj. (A.6) 



We only consider 5 > 2. Note that we have used the definition of the reduced Planck mass, Eq. (2.7). 
The term in the square brackets equals e/Arri^^^ (as defined in Sec. 2). 



4-fermion operators for fx -^ see and fj, —^ e conversion in nuclei 

In addition to the magnetic moment type operator Eq. (A. 2), the following 4-fermion operator contributes 
to the fi -^ eee process: 



J 



L 



V 



e7„PL/i 67" {{A^ + glAz + Ab) Pl + (^7 + 9rAz) Pr] e + h.c. 



(A.7) 



Here g^ff\ = PIl(r) ~ QfS>vii 9w (ffi = —1/2 + sni 6w and g^^ = sin^iy)- The coefficients A^, Az, 
and Ab correspond to contributions from photon penguin, Z-penguin, and box diagrams, respectively. 
Explicitly, 



A. 



'327r2 



771 



w 



Ey^HV'HF.imr/ml.), 



(A.8) 



n.,J 



647r2 



771 



W 



Y^n^llFzirnT Im 



W) 



n,J 



+ E ylHVrvrX:^Gzimr/ml,,m%^Vml.)\, 



,0-)2, 



n,m,j 



(A.9) 
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Ab 

F^{x) 

Fz{x) 

Gz{x,y) 

GB{x,y) 



647r2 



m 



W 



E yenVllVl^VilGBimt^'/m^w, m'^'/mlr), 



n,m.,i,] 



X {12 - llx - Sx^ + Tx^ + 2x{x^ - lOx + 12) In x} 
12(1 -x)4 ' 

5x(l — X + xlnx) 
(l-x)2 ' 

1 Jx2(l — y)lnx y^(l — x)\ny\ 



x-y 

1 
x-y 

-2xy 



1-x 



xy 
1 + — 



1 



+ 



1-y 

x^ Inx 



(A.IO) 

(A.ll) 
(A.12) 
(A.13) 



1 



1 xlnx 

+ 



1 — X ' (1 — x)2 1 — 2/ (1 — y)^ 

1 ylny 



1 — X (1 — x)2 1 — y (1 — y)^ 
The branching ratio for /i ^> eee process is 



(A.14) 



BR(^ -^ eee) 



1 



8GI 



\A^ + g'JiAzl^ + 2\A^+ glAz + Ab\^ + 32 {eA^l^ In 



mt. 11 



m" 



+4e^H {3^; + {9% + 2gi)A^ + 2A*s} + h.c. 



(A.15) 



In addition to the magnetic moment type operator Eq. (A. 2), the following 4-fermion operator contributes 
to fi ^ e conversion in nuclei: 



q=u,d V / 



B 



B 



1287r2 mfy 



^ E K^^M^ |K,.|^G^(m(^')V-^v.,</-'v.), 



J^.'j 



B% = 

The /Li ^ e conversion rate is 
R(/i ^ e) = 



9' 1 



327r2 



m 



W rJ 



EK^^^^^GbK^^V-^vf-O). 



(A.16) 
(A.17) 
(A.18) 



"J 



a3 z,4jF(-m2)|2„, 



47r2 Z 



^— capt 



{2Z + N)(^^4^Az + B^ 



+ {Z + 2N) m±^Az + bU - Z{2eAR + A^) 



(A.19) 



where r^_capt is the muon capture rate in nuclei of interest [46], Z and N are the proton and neutron 
numbers, respectively, F{—mV) is the nuclear form factor and Zgfj is the nuclear effective charge [47]. 
Numerically, these nuclear parameters are r^_capt = 1.7x10"^^ GeV (4.6x10"^^ GeV), Z^s = 17.6 (11.6), 
and \F{-ml)\ = 0.535 (0.64) for ||Ti (f^Al) [46, 47]. 

Summing over the KK states following the steps previously outlined, we obtain 



A 



Az 



e'5' 1 



11. —e^u r/*. 

647r2m2^ ^^'^^^^ 



m 



vj 



A^T-atm |6 ' 3 



7 1 A A2 
- + - In ^ 



m 



W 



5-2 



(A.20) 
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647r2 
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Eueju;^ 



m 



^3 



^"T-atm ' 



(A.21) 



(A.22) 

(A.23) 
(A.24) 
(A.25) 



where Js = (1 — 6/2)"^ /q /q {zw) '"^ ^ \\i{z/w)/{z — w) dz dw, which is the number of order 1. Numerically, 
Js = 0.23, J4 = 0.43, J5 = 0.55, and Jq = 0.63. 

Note that, unlike the /i — > 67 case, most of the amplitudes depend not only on the neutrino mass- 
squared difference (which is directly measured by neutrino oscillation experiments), but also on the 
magnitude of the neutrino mass-squared. 



Magnetic moment operator for muon g — 2 

The anomalous magnetic moment ((5a^) of the muon is defined as the coefficient of the effective operator 



C 



Am, 



-5a^ fiaa/BfJ-F 



a/3 



(A.26) 



The contribution to muon anomalous magnetic moment induced by massive KK modes is given by 



6a„ 



fix) 



g^ m^ 



167r2 






6/' 



10 - 43x + 78x2 _ 4g^3 _^ 4^4 _^ ]^g^3 1^ ^ 
12(1 -3;)4 ■ 



(A.27) 
(A.28) 



Note that we have subtracted the SM {W, i^)-loop in order to calculate the new physics contribution. It 
is useful to separate the sum into the "massless" part (we neglect the effect of the small active neutrino 
masses) and the "massive" part. Using the unitarity of V^^^.j and the fact that /(O) = 5/6, we can rewrite 



5a,, 



5f2 m? 



167r2 



m 



W 



fEKA G,{mr/ml,) 



(A.29) 



nj 



Here, the function G^{x) = f{x) — 5/6 is defined by Eq. (A. 3). Summing over the KK states (see 
Eq. (A.6)) we obtain Eq. (4.17). 
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